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We study the effect of a potential fourth fermion generation on the upper and lower Higgs boson 
mass bounds. This investigation is based on the numerical evaluation of a chirally invariant lattice 
Higgs- Yukawa model emulating the same Higgs-fermion coupling structure as in the Higgs sector 
of the electroweak Standard Model. In particular, the considered model obeys a Ginsparg-Wilson 
version of the underlying SU(2)^ x U(l)y symmetry, being a global symmetry here due to the 
neglection of gauge fields in this model. We present our results on the modification of the upper 
and lower Higgs boson mass bounds induced by the presence of a hypothetical very heavy fourth 
quark doublet. Finally, we compare these findings to the standard scenario of three fermion 
generations. 
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1. Introduction 

The Sakharov explanation for the matter anti-matter asymmetry of the universe suffers from 
the CP-violating phase of the Standard Model (SM3) falling short by at least 10 orders of mag- 
nitude. In addition to this concern the Sakharov picture demands a first order electroweak phase 
transition, which is also objected in the framework of the SM3. However, both of the above caveats 
might be addressable [jjj 0] by the inclusion of a new fourth fermion generation into an extended 
version of the Standard Model (SM4). Despite the arguments against the existence of a fourth 
fermion generation such a scenario nevertheless remains attractive for two reasons. Firstly, there is 
a strong conceptual interest, since a new fermion generation would need to be very heavy, leading 
to rather large Yukawa coupling constants and thus to potentially strong interactions with the scalar 
sector of the theory. Secondly, it has been argued ^ (and the references therein) that the fourth 
fermion generation is actually not excluded by electroweak precision measurements, thus leaving 
the potential existence of a new fermion generation an open question. 

In our contribution, however, we do not present any statement arguing in favour or disfavour 
of a new fermion generation. Here, we simply assume its existence and focus on the arising conse- 
quences on the Higgs boson mass spectrum. With the advent of the LHC this question will become 
of great phenomenological interest, since the Higgs boson mass bounds, in particular the lower 
bound, depend significantly on the heaviest fermion mass. Demonstrating this effect will be the 
main objective of the present work. 

Due to the large Yukawa coupling constants of the fourth fermion generation a non-perturbative 
computation is highly desirable. For this purpose we employ a lattice approach to investigate the 
strong Higgs-fermion interaction. In fact, we follow here the same lattice strategy that has already 
been used in Ref. [||] for the non-perturbative determination of the upper and lower Higgs boson 
mass bounds in the SM3. This latter approach has the great advantage over the preceding lattice 
studies of Higgs-Yukawa models that it is the first being based on a consistent formulation of an 
exact lattice chiral symmetry [Q, which allows to emulate the chiral character of the Higgs-fermion 
coupling structure of the Standard Model on the lattice in a conceptually fully controlled manner. 

2. Numerical Results 

In order to evaluate the Higgs boson mass bounds we have implemented a lattice model of 
the pure Higgs-fermion sector of the Standard Model. To be more precise, the Lagrangian of the 
targeted Euclidean continuum model we have in mind is given as 

- 1 1 2 - 

L HY =f?t' + b'?b' + -dn(p*d IJt (p + -ml(p*(p + X (<p» +yy [f ,b') L (pb' R +y t > (f,b') L (pt R 

+ c.c. of Yukawa interactions, (2. 1) 

where we have constrained ourselves to the consideration of the heaviest quark doublet, i.e. the 
fourth generation doublet, which is labeled here (t',b'). This restriction is reasonable, since the 
dynamics of the complex scalar doublet q> ((p = z'T2<p*, T; : Pauli-matrices) is dominated by the 
coupling to the heaviest fermions. For the same reason we also neglect any gauge fields in this 
approach. The quark fields nevertheless have a colour index which actually leads to N c = 3 identical 
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copies of the fermion doublet appearing in the model. However, for a first exploratory study of the 
fermionic influence on the Higgs boson mass bounds we have set 7V C to 1 for simplicity. 

The actual lattice implementation of the continuum model in Eq. ( |Q| ) has been discussed in 
detail in Ref. []|]. Since the Yukawa interaction has a chiral structure, it is important to establish 
chiral symmetry also in the lattice approach. This has been a long-standing obstacle, which was 
finally found to be circumventable by constructing the lattice equivalent of <j? as well as the left- 
and right-handed components of the quark fields t' LR , b' LR on the basis of the Neuberger overlap 
operator [Q BJ]. Following the proposition in Ref. [Q] we have constructed a lattice Higgs- Yukawa 
model with a global SU (2)l x U (l)y symmetry. 

The fields considered in this model are the aforementioned doublet cp as well as N c quark dou- 
blets represented by eight-component spinors = i = \,...,N C . With £^( ov ) denoting 
the Neuberger overlap operator the fermionic action S> can be written as 

N c 

Sf = £ = @ {ov) +P+^ tiag(yt>,y v )P++P- diag(^,y y )^_, (2.2) 

i=l 

where y t ip denote the Yukawa coupling constants and the scalar field cp x has been rewritten here as 
a quaternionic, 2x2 matrix (j)£ = (<p x , (p x ), with x denoting the site index of the Lr s x L r lattice. The 
left- and right-handed projection operators P± and the modified projectors P± are given as 

P ± = ^p, P± = ^, 7 5 = 75 (l-^ (ov) ), (2.3) 

with p being the radius of the circle of eigenvalues in the complex plane of the free Neuberger 
overlap operator [BJ. 

This action now obeys an exact global S\J(2)l x U(l)y lattice chiral symmetry. For Q.^ G 
SU(2) and £6R the action is invariant under the transformation 

Y -> U Y P + y + U Y 0. L P- V, W -> W p + a l U Y + W- u y > ( 2 - 4 ) 
-»• UrQ&l, 1 " -> £Il^U~y (2-5) 

with the compact notation Uy = exp(z£F) denoting the respective representations of the global 
hypercharge symmetry group U(l)y for the respective field it is acting on. In the continuum limit 
Eq. (^4]-^5]) eventually recover the (here global) continuum SU(2)i x U(l)y chiral symmetry. 

Finally, the purely bosonic part Sy of the total lattice action S = Sp + S v is given by the usual 
lattice <p 4 -action 



S f = £ ^W x ^i<Px + ~m 2 o<Px<Px + A (<pl<p x f , (2.6) 



with the bare mass mo, the forward difference operator in direction jj., and the bare quartic 



coupling constant A . For the practical lattice implementation, however, a reformulation of Eq. ( pLq ) 
in terms of the hopping parameter K and the lattice quartic coupling constant X proves to be more 
convenient. It reads 

5* = [<t>x+» + <£*- M ] + £ ®t* x + X £ -N c f, (2.7) 

x,n x x 
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and is fully equivalent to Eq. (2.6). This alternative formulation opens the possibility of explicitly 
studying the limit A = °° on the lattice. The aforementioned connection can be established through 
a rescaling of the scalar field & x € 1R 4 and the involved coupling constants according to 



'2k 



® 2 x + i& x 



1-2AU-8K 

K 



(2.8) 



Due to the triviality of the Higgs sector the targeted Higgs boson mass bounds actually depend 
on the non-removable, intrinsic cutoff parameter A of the considered Higgs-Yukawa theory, which 
can be defined as the inverse lattice spacing, i.e. A = I /a. To determine these cutoff dependent 
bounds for a given phenomenological scenario, i.e. for given hypothetical masses of the fourth 
fermion generation, the strategy is to evaluate the maximal interval of Higgs boson masses attain- 
able within the framework of the considered Higgs-Yukawa model being in consistency with this 
phenomenological setup. The free parameters of the model, being the bare scalar mass mo, the 
bare quartic coupling constant A, and the Yukawa coupling constants y t i y thus have to be tuned 
accordingly. The idea is to use the phenomenological knowledge of the renormalized vacuum ex- 
pectation value v r /a = 246 GeV of the scalar field (p as well as the hypothetical fourth generation 
quark masses m t i % y to fix the bare model parameters for a given cutoff A. 

In lack of an additional matching condition a one-dimensional freedom is left open here, which 
can be parametrized in terms of the quartic coupling constant A. This freedom finally leads to the 
emergence of upper and lower bounds on the Higgs boson mass. As expected from perturbation 
theory, one also finds numerically ^ that the lightest and heaviest Higgs boson masses are obtained 
at vanishing and infinite bare quartic coupling constant, respectively. The lower mass bound will 
therefore be obtained at X = 0, while X = oo will be adjusted to derive the upper bound. 

Concerning the hypothetical masses of the fourth fermion generation quarks, we target here a 
mass degenerate scenario with m t >/a = my /a = 700GeV, which is somewhat above its tree-level 
unitarity upper bound [g]. However, we are currently also investigating a set of other mass settings 
to study in particular the quark mass dependence of the Higgs boson mass bounds. 
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Table 1: The model parameters underlying the lattice calculations performed in this study are presented. 
The setting A = (A = °°) is employed for deriving the lower (upper) Higgs boson mass bound. Depending 
on the lattice volume the available statistics ranges from Nconf = 1,000 to Nconf — 20,000. 
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For the eventual determination of the cutoff dependent Higgs boson mass bounds several series 
of Monte-Carlo calculations have been performed at different values of A and on different lattice 
volumes as summarized in Tab. [[[ In order to tame finite volume effects as well as cutoff effects 
to an acceptable level, we have demanded as a minimal requirement that all particle masses m = 
mji,m t i,my in lattice units fulfill m < 0.5 and m-L st > 3.5, at least on the largest investigated lattice 
volumes. Assuming the Higgs boson mass m H to be around 500 — 750GeV this allows to reach 
cutoff scales between 1500GeV and 3500GeV on a 24 3 x 32-lattice. However, despite this setting 
strong finite volume effects are nevertheless expected induced by the massless Goldstone modes. It 
is known that these finite size effects are proportional to 1 /L 2 S at leading order. An infinite volume 
extrapolation of the lattice data is therefore mandatory. 




(d) 




0.004 0.006 
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Figure 1: The finite volume data of the renormalized vacuum expectation value v r (a,d), the Higgs boson 
mass m,H (b,e), and the degenerate quark mass m,i — my (c,f), as obtained from the lattice calculations 
specified in Tab. [j], are plotted versus 1/1,2. Xhe 

upper (lower) row corresponds to the setting A = (A = °°). 
The infinite volume extrapolation is performed by fitting the data to a linear function. Due to the observed 
curvature arising from the non-leading finite volume corrections only those data with L s > 16 have been 
respected by the linear fit procedures. 



The finite volume data of the renormalized vacuum expectation value v r , the Higgs boson 
mass ma, and the degenerate quark mass m t i = my resulting from the calculations specified in 
Tab. [j] are presented in Fig. []]. For the details about the determination of the latter observables the 
interested reader is referred to Ref. [Qj. Here it is only stated that the renormalization constant Zq 
entering the renormalization of the scalar field q> has been derived from the Goldstone propagator, 
the Higgs boson mass mn has been taken from the Higgs propagator, and the quark masses have 
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been computed from the time correlation function of the respective fermionic fields. These lattice 
data are plotted versus \/L 2 s and extrapolated to the infinite volume limit by means of a linear fit 
ansatz according to the aforementioned leading order behaviour. Due to the observed curvature 
arising from the non-leading finite volume corrections only those data with L s > 16 have been 
respected by the linear fit procedures. One finds that the intended infinite volume extrapolation can 
indeed reliably be performed thanks to the multitude of investigated lattice volumes reaching from 
12 3 x 32 to 24 3 x 32-lattices here. 

The quality of the tuning procedure intending to hold the quark masses constant can then 
be examined in Fig. |2|b displaying the results of the infinite volume extrapolation of m t i. In the 
considered SM4 scenario the fluctuation of the quark mass has been constrained to roughly m t i = 
my = 676 ±22 GeV. For later comparisons with the corresponding SM3 scenario we also present 
the analogous summary plot of our earlier investigations [^] of the latter setup where the degenerate 
quark masses have been fixed to approximately m t = mj, = 173 ±3 GeV as demonstrated in Fig. 0a. 
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Figure 2: The infinite volume extrapolations of the degenerate quark masses observed in the lattice calcula- 
tions specified in Tab. |]are presented versus the cutoff parameter A. In the SM3 scenario (a) the fluctuation 
of the quark mass has been constrained to m t = mt = 173 ± 3GeV, while m t i — my — 676 ± 22 GeV is 
adjusted in the SM4 scenario (b). 



The infinite volume results of the Higgs boson masses are finally presented in Fig. ||b. The 
numerical data for the upper mass bound have moreover been fitted with the analytically expected 
functional form of the cutoff dependence of the upper Higgs boson mass bound derived in Ref. [Q]. 
It is given as 



up 

^ =A m .[log(A 2 /M 2 )+B» 
a 



-1/2 



(2.9) 



with A m , B m denoting the free fit parameters and jit being an arbitrary scale, which we have cho- 
sen as jU = 1 TeV here. One learns from this presentation that the obtained results are indeed in 
good agreement with the expected logarithmic decline of the upper Higgs boson mass bound with 
increasing cutoff parameter A. 
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The reader may want to compare these findings to the upper and lower Higgs boson mass 
bounds previously derived in the SM3. The lattice results corresponding to that setup have been 
determined in Ref. [||] and are summarized in Fig. ||a. The main finding is that especially the lower 
bound is drastically shifted towards larger values in the presence of the assumed mass-degenerate 
fourth quark doublet. From this analysis it can be concluded that the usually expected light Higgs 
boson seems to be incompatible with a very heavy fourth fermion generation. 
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Figure 3: Upper and lower Higgs boson mass bounds are shown foiN c = 1, m t = m\ = 173 ± 3GeV (a) and 



N c = l, m t i = my = 676 ± 22 GeV (b). Both upper bounds are each fitted with Eq. ( 2.9 ). The lower bound 
in (a) is also compared to a direct analytical computation depicted by the solid line as discussed in Ref. [||. 
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